See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/299282957 
Automaty komorkowe jako modele wzrostu populacji 


Chapter - January 2007 


CITATIONS READS 
0 201 
1 author: 


Fo WOJCIECH BORKOWSKI 
University of Warsaw 
93 PUBLICATIONS 696 CITATIONS 


SEE PROFILE 


Some of the authors of this publication are also working on these related projects: 


Project Languages & linguistics View project 


Project Open Science & Creative activity View project 


All content following this page was uploaded by WOJCIECH BORKOWSKI on 05 July 2018. 


The user has requested enhancement of the downloaded file. 


ResearchGate 


Wojciech Borkowski 


Automaty komórkowe jako modele wzrostu populacji 

(In book: Modelowanie matematyczne i symulacje komputerowe w naukach 
społecznych Edition: Warszawa, Publisher: ACADEMICA 
https://www.researchgate.net/publication/299282957_Automaty_komorkowe_jako_mo 


dele_wzrostu_populacji ) 


STRESZCZENIE: W tradycji matematycznej systemy dynamiczne modeluje się zwykle „„makroskopowo” 
uwzględniając jedynie najważniejsze globalne właściwości i wprowadzając wiele upraszczających założeń. 
Niestety jednym z najpopularniejszych jest pominięcie zróżnicowania przestrzennego i znaczenia lokalnych 
zmian, które zarówno dla wielu dziedzin biologii jak i dla nauk społecznych są jednym z ważniejszych 
aspektów. Automaty komórkowe są najprostszym sposobem alternatywnego modelowania, „„mikroskopowo” 
uwzględniającym zróżnicowanie przestrzenne w systemach. Trywialny z matematycznego punktu widzenia 
proces zasiedlanie obszaru przez populację jest dobrym pretekstem by zdefiniować zasady i pokazać efekty 
działania kilku automatów komórkowych, które odwołując się do przykładów z rzeczywistości biologicznej i 
społecznej, symulują to zjawisko z rosnącą dokładnością, aż po osiągnięcie jednoczesnej zgodności z modelem 
„makroskopowym”. Wszystkie opisane modele należą do jednej ,,rodziny”, której szczególnym przypadkiem jest 
najbardziej znany automat komórkowy: „Gra w Życie” („Game of Life”) — sztandarowy przykład systemu, który 
mimo prostoty reguł mikroskopowych przejawia niesłychane bogactwo emergentnych zachowań 


makroskopowych. 


W klasycznych modelach matematycznych układów dynamicznych przyjmuje się 
różne upraszczajace założenia umożliwiające analizę modelu metodami analizy 
matematycznej. Niestety jednym z głównych założeń jest uznanie, że przestrzeń nie ma 
znaczenia — wszelkie lokalne różnice przestrzenne badanego układu „uśredniają się” i mogą 
być zignorowane. Zapewne w wielu przypadkach jest to słuszna intuicja, ale łatwo wskazać 


sytuacje gdy może ona prowadzić do błędnych wniosków. 


Weźmy za przykład ogólny model wzrostu dowolnej populacji, analogiczny do klasycznego 
równania logistycznego, w którym nowych indywiduów przybywa proporcjonalnie do liczby 
już istniejących, ale pod warunkiem, że populacja nie jest zbyt liczna. Im więcej indywiduów 
już mamy, tym wolniej ich przybywa, aż do momentu, gdy liczba nowo powstałych, czy 


przybyłych zaledwie równoważy śmiertelność lub emigracje. 


Zauważmy, że opis ten jest wybitnie makroskopowy - traktuje populacje jako całość nie 
zastanawiając się nad rozmnażaniem i śmiercią konkretnych osobników. Co więcej cała 
populacja jest wyrażana jedynie jako liczba w zakresie 0 do 1. Wiemy, w jakim stopniu 
populacja wypełnia swoje środowisko — w 5, czy w 50 czy w 99%, ale nie wiemy, jaką 
konkretną liczbę osobników to oznacza. Nie wiemy też czy osobniki są rozrzucone 
równomiernie, parami, czy może tworzą jedno wielkie stado. Nie wiemy czy możemy zgodzić 
się na populację zapełniającą środowisko w 0,000001%? Jeśli będzie to populacja dużych 
zwierząt to zapewne nie - ćwierć słonia czy pół hipopotama to za mało by populacja mogła 
przeżyć! Nawet jeśli liczba ta oznacza kilka osobników to nie wiemy jak są rozmieszczone i 


czy będą w stanie kiedykolwiek się spotkać by wydać na świat potomstwo... 


Czy możemy modelować układy dynamiczne, dla których istotne są lokalne oddziaływania 
i zróżnicowanie przestrzenne? Oczywiście jest to możliwe, ale potrzebny jest model zupełnie 
innego typu. Musimy porzucić spojrzenie generalizujące badany układ i pochylić się nad 
indywidualnymi losami budujących go obiektów, czyli wyjść od „mikroskopowego” 
spojrzenia na zjawisko. Może to oznaczać konieczność uwzględnienia pojedynczych 


osobników populacji i ich lokalnych oddziaływań. 


Najprostszy automat komorkowy 


W przypadku modelu wzrostu populacji, czy raczej „zasiedlania obszaru”, będą nas 
interesować oddziaływania prowadzące do przeżywania, rozmnażania się i śmierci 
osobników. Jak w każdym modelu musimy zgodzić się na pewne uproszczenia. Na początek 
przyjmiemy, że interesująca nas przestrzeń jest jednowymiarowa. Podzielimy ją na odcinki 
równej długości, z których każdy może być „zamieszkany? przez najwyżej jedno 
indywiduum. Innymi słowy mamy ciąg obiektów - komórek, które mogą być zajęte lub nie. 
Możemy sobie wyobrazić, że interesuje nas populacja ptaków gniazdujących nad potokiem. 
Tak jak to często bywa w przypadku ptaków terytorialnych ważne są tylko gniazdujące pary. 
Każdy samiec zajmuje swój odcinek potoku i nie dopuszcza do niego innych. Dzięki temu 
może zwabić do swego gniazda jedną z wędrujących po okolicy samic i wychować pisklęta. 
W następnym roku synowie zajmują sąsiednie terytoria, o ile są wolne, albo podobnie jak 
córki odlatują „w siną dal” w poszukiwaniu innego potoku i przestają nas interesować. 
Każdego roku coraz większa liczba par będzie gniazdować w koloni rozciągniętej wzdłuż 


potoku. 
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Rysunek 1: Jeden z najprostszych jednowymiarowych automatów komórkowych jako model terytorialnej 


populacji ptaków. W wersji dosłownej oraz przedstawiony nieco bardziej abstrakcyjnie. 


Jeśli puste terytoria oznaczymy zerami, a zajęte jedynkami to otrzymamy jeden z 
najprostszych automatów komórkowych (rys. 1). Posiada on wszystkie niezbędne 


właściwości tej klasy modeli: 


e Przestrzeń jest podzielona równomiernie na regularne komórki mogące przejmować 
co najmniej dwa różne stany. Tutaj mamy nadrzeczne terytoria, które mogą być 


zasiedlone (komórki z 1) lub nie (komórki z 0). 


e Stany poszczególnych komórek zmieniają się zgodnie z lokalnymi regułami, czyli w 
zależności od stanu danej komórki i jej sąsiadów. W przykładzie reguła zmiany stanu 
brzmi „komórka zostaje zasiedlona, jeśli sąsiaduje z zasiedloną komórką”, jest zatem 


niezwykle prosta. 


e Czas płynie dyskretnie - stan całego automatu, czyli zbiór stanów wszystkich jego 
komórek zmienia się krokami. Dla naszej populacji ptaków jeden krok to jeden sezon 


lęgowy czyli jeden rok. 


Zasiedlanie kempingu nad rzeką 


Zajmijmy się teraz odrobinę trudniejszym przykładem. Wyobraźmy sobie, że zamiast 
nadrzecznych zarośli mamy rozłożony wzdłuż wybrzeża kemping podzielony na stanowiska 
równej długości — to będą nasze komórki. Na każdym stanowisku może stać najwyżej jeden 
namiot. Przyjmijmy, że osoby wybierające miejsce na swój namiot nie chcą by stał samotnie, 
bo wtedy wyróżnia się i jest bardziej narażony na splądrowanie przez grasujących 
„rabusiów ”, nie chcą też by namiot stał całkowicie otoczony innymi, bo w tłoku nikt nie czuje 
się komfortowo. Przy wyborze miejsca biwakowania kierują się zatem następującymi 
regułami: 
e Jeśli chcesz postawić namiot to wybierz miejsce, które sąsiaduje dokładnie z jednym 
już stojącym namiotem 
e Jeśli już biwakujesz to postanawiasz pozostać na następny dzień jeśli twój namiot ani 
nie stoi samotnie ani nie jest otoczony z obu stron. 
Reguły te można uprościć i powiedzieć, że „komórka może zostać lub pozostać zasiedlona 
jeśli sąsiaduje z dokładnie jedną zasiedloną komórką”. 
Przyjmiemy też, że recepcja jest czynna tylko godzinę dziennie. Otwierana jest godzinę po 
przyjeździe jedynego autobusu kursującego do tej miejscowości i zamykana godzinę przed 


jego odjazdem z powrotem do miasta. Zatem decyzję o zameldowaniu lub wymeldowaniu, 


albo o zmianie miejsca biwakowania trzeba podjąć na podstawie stanu sprzed otwarcia 
recepcji, a potem ustawić się w kolejce. Choć w dzisiejszych czasach to nieco sztuczne 
założenie', potrzebujemy go, by móc nadal myśleć, że czas w modelu płynie dyskretnymi 


skokami, a nie w sposób ciągły. 
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Rysunek 2: Jednowymiarowy automat komórkowy jako model zasiedlania nadrzecznego kempingu. 
Namioty i jedynki pokazują komórki zasiedlone, osoby z podniesionymi rękami pokazują komórki zdatne 
do zasiedlenia w następnym kroku — dniu. 

Możemy prześledzić historię kempingu na kartce papieru, ale lepiej za pomocą komputera. 
(rys. 2). Nie bez powodu tego typu modele, choć wymyślone już w latach 40-tych XX wieku 
przez Stanisława Ulama i Johna von Neumana, zaczęły być powszechniej badane dopiero od 
początku lat 70-tych, gdy komputery stały się na tyle tanie by rozpowszechnić się w 
środowisku naukowym. Zazwyczaj badanie automatów komórkowych wymaga żmudnego 
powtarzania ogromnej liczby prostych obliczeń. Ludzie słabo sprawdzają się w takich 


zadaniach, za to komputery znakomicie! 


Jeśli rozpoczniemy od pojedynczego namiotu, to historia stanów automatu rysuje regularną 
figurę — coś w rodzaju powycinanego trójkąta”. Wynik ten jest nieco zaskakujący - trudno by 
było przewidzieć takie zachowanie tylko na podstawie znajomości reguł. Jeszcze trudniej 
przewidzieć jak będą wyglądały kolejne stany automatu, jeśli wystartuje on z jakiejś 
konkretnej konfiguracji jedynek — namiotów. Dzięki implementacji modelu w komputerze 
możemy jednak prześledzić wiele różnych konfiguracji i spróbować wyciągnąć jakieś ogólne 


wnioski. W tym raczej prostym przypadku wystarczy sprytnie wypełniony arkusz 


' Ale całkiem realistyczne jeszcze w latach 80-tych, gdy biwakująca młodzież przemieszczała się niemal 
wyłącznie pociągami i PKS-ami, a praca nie była cennym dobrem, lecz „dopustem”. 
> Nawiasem mówiąc jest to przybliżenie tzw. trójkąta Sierpińskiego — figury fraktalnej odkrytej przez 


matematyków w zupełnie innych okolicznościach. 


kalkulacyjny, ale zwykle do badania automatów komórkowych używa się specjalnie do tego 
celu napisanych programów”. Do wygenerowania konfiguracji początkowych posłużymy sie 
losowaniem, a ściślej mówiąc tzw. generatorem liczb pseudolosowych. Dzięki niemu 
będziemy mogli ze ściśle określonym prawdopodobieństwem wylosować stan początkowy dla 
każdej komórki, czyli zadecydować czy w komórce będzie 0 czy 1. Prawdopodobieństwo 
równe % albo 50% oznacza, że jedynki i zera są równie prawdopodobne, a 
prawdopodobienstwo 1⁄4 czyli 25% oznacza, że jedynki i zera będą wylosowane w proporcji 
1:3. Żeby uniknąć kilku problemów, o których nie pora teraz mówić, będziemy ze znacznie 
większego obszaru wylosowanego z danym prawdopodobieństwem obserwować fragment o 


tej samej długości (50) i dokładnie tej samej początkowej liczbie jedynek. 


> Wiele z nich jest łatwo i za darmo dostępne w Internecie. Np. znakomity program MCell autorstwa Mirosława 


Wójtowicza można ściągnąć ze strony http://www.mirekw.com/ca/ 
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fragmenty trójkąta Sierpińskiego pojawiające się wtedy, gdy w obrębie większego obszaru, 


który został wyzerowany, pozostanie samotna jedynka. 


Wydaje się też, że zmniejszona lub zwiększona liczba jedynek na początku wpływa najwyżej 
na kilka pierwszych kroków, dalsze są już nieodróżnialne — mogłyby pochodzić z dowolnej 
historii. Zatem automat dąży do tego by jego wypełnienie jedynkami było zamknięte w 
pewnych granicach. Możemy przyjrzeć się temu dokładniej. Wystarczy policzyć liczbę 
jedynek w każdym kroku modelu. Otrzymane dane umożliwiają ilościowe porównanie wielu 
różnych historii startujących z różnych konfiguracji (rys. 4) 1 stwierdzenie, że chociaż liczba 
jedynek zmienia się nieregularnie (oscyluje) to faktycznie rzadko wykracza poza środkową 
część zakresu (tu 15-35 jedynek), przynajmniej dopóty, dopóki nie zaczynamy od 


„ekstremalnie” małej lub dużej liczby jedynek. 
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Rysunek 4: Zmiana liczby jedynek (N) w 50 dniach (krokach) w 50 komórkowych wycinkach automatu 
modelujacego „kemping”. Górny wykres prezentuje przebiegi dla historii a, b, c z rysunku 3. Dolny 
wykres prezentuje po 5 przebiegów analogicznie startujących z losowych konfiguracji 13, 25 i 37 


komórek. 


Dwuwymiarowy, komorkowy model wzrostu 


Co stanie się, gdy nasz model kempingu zechcemy rozszerzyć na dwa wymiary, czyli zająć 
się np. wielkim polem namiotowym przed koncertem rockowym „Przystanek Woodstock” 
albo w Jarocinie, albo może wielkim spotkaniem młodzieży katolickiej w Lednicy? W 
najprostszym przypadku komórki z odcinków staną się kwadratami. Przestrzeń można by też 
podzielić na trójkąty równoboczne lub sześciokąty, ale łatwiej myśleć o modelu, w którym 
przestrzeń komórek będzie przypominać plansze do gry w warcaby — chociaż może znacznie 
większą, a teoretycznie nawet nieskończona. Oczywiście każda komórka będzie teraz miała 
maksymalnie nie dwu bezpośrednich sąsiadów, ale nawet ośmiu (rys. 5a). Reguła zmiany 
stanu z automatu jednowymiarowego może zostać zachowana, lecz przestaje być jedyną 
możliwą regułą implementującą zachowanie typu „nie za dużo, nie za mało”. Jest za to 
najbardziej skrajną! Kogoś, kto w takiej sytuacji toleruje posiadanie tylko jednego sąsiada 


uznalibyśmy raczej za samotnika. 


Powinniśmy jednak zacząć od sprawdzenia działanie takiej najprostszej reguły. Na rysunku 
(rys. 5b) widać kilka kolejnych kroków automatu zainicjowanego jedną komórką. Jeśli 
wyobrazimy sobie, że kolejne kroki leżą jak kolejne warstwy jedna pod drugą, to widać 
wyraźną analogię z wersją jednowymiarową. Gdy jednak zastosujemy inicjację losową to 
porządek znika. Otrzymujemy obraz (rys. 5c) pozbawiony niemal rozpoznawalnych wzorców, 
a do tego bardzo zmienny w czasie. Na wykresie zmian liczby ,.namiotow” (rys. 5d) widać, że 
początkowo bardzo znaczne wahania z czasem gasną i rozmiar populacji stabilizuje się. 
Jednak już dokładna obserwacja rysunku pozwala się zorientować, że mało która z komórek 
może pozostać zasiedlona do następnego kroku. Większość „namiotów” w każdym kroku 
stanowią te nowo postawione. Przy tak małej tolerancji na tłok jedynie wielka potencjalna 
„płodność” komórek (aż do 8 z 1) pozwala na utrzymanie zasiedlenia. Nie jest to więc zbyt 
realistyczny model, ale na szczęście możemy wypróbować także inne reguły optymalizujące 


zasiedlanie. 


Choćby taką, że wypełniona komórka (,,zywa’’, reprezentowana przez 1) może pojawiać się i 
pozostawać w miejscu gdzie sąsiaduje z dwiema lub trzema zajętymi komórkami. To 
założenie może nieco arbitralne dla pola namiotowego, dla modelu biologicznego, od którego 
zaczęliśmy byłoby całkiem rozsądne. W końcu większość organizmów, które da się zobaczyć 
gołym okiem potrzebuje dwojga rodziców i jest w stanie tolerować obecność rodzeństwa. Czy 
znając działanie poprzednich modeli możemy pokusić się o jakieś przewidywanie 


zachowania? Zatrzymajmy się na chwilę, i pomyślmy, co można uzyskać zasiedlając model 


parą „pierwszych rodziców”, a co gdy zasiedlenie będzie losowe np. z prawdopodobieństwem 


zasiedlenia komórki wynoszącym 5%? 
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Rysunek 6: Automat komórkowe o regule „żyj gdy masz 2 lub 3 żywych sąsiadów”. A) Dwa kroki dla 
pary „żywych” komórek rozmieszczonych po przekątnej - tworzącej tzw. oscylator oraz pierwsze 
przekształcenie pary ustawionej w jednej linii — rozpoczynającej wzrost rozety. B) Ekspandująca rozeta 
po stu kilkudziesięciu krokach. C) Ten sam fragment po kolejnych kilkuset krokach. D) Wiele rozetek 
krótko po inicjacji losowej z prawdopodobieństwem 1%. E) ...i po stu kilkudziesięciu krokach 

Tym razem model zachowuje się dosyć rozsądnie”. Przy inicjacji parą rodziców 
otrzymujemy ekspandującą rozetę, która po odpowiednio długim czasie wypełnia całą 
dostępną przestrzeń skomplikowanym i dynamicznym wzorem, ale z zachowaniem 
początkowej symetrii (rys. 6bc). Gdy populacja wystartuje od losowego rozmieszczenia 
zasiedlonych komórek, te pojedyncze natychmiast giną, pary sąsiadujące pionowo i poziomo, 
a także znacznie rzadsze konfiguracje 3 i 4 sąsiadujących komórek dają początek rozetom. 


Natomiast pary sąsiadujące po przekątnej giną, lecz tworzą nowe pary komórek sąsiadujące 
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po przekątnej (rys. 6a), a w kolejnym kroku są ponownie odtwarzane w pierwotnym miejscu i 
tak „w nieskończoność”. Takie lokalne, dynamiczne konfiguracje odtwarzające stan 
początkowy co 2 lub więcej kroków zwykło nazywać się „oscylatorami”. W kolejnych 
krokach obszary ekspansji poszczególnych rozet prędko się zazębiają, pochłaniają pary 
oscylujące i resztę dostępnej przestrzeni. Po pewnym czasie otrzymujemy „kaszkę” bez 
żadnych rozpoznawalnych wzorców (rys. 6cd), która co prawda wciąż zmienia się lokalnie, 


ale pozostaje w stanie równowagi o niezmiennej gęstości zasiedlonych komórek. 


Analogicznie zachowają się modele, w których zwiększymy tolerancję na tłok do 4 czy 
więcej „żywych” komórek w sąsiedztwie. Jedyną różnicą będzie gęstość komórek we wnętrzu 
rozet, i gęstość zasiedlenia w stanie równowagi. W skrajnym przypadku, gdy komórki tolerują 
obecność nawet 8 „żywych” sąsiadów przestrzeń zostanie oczywiście wypełniona całkowicie 


i dynamika modelu ustanie. 


Realistyczny model zasiedlenia 


Model dwuwymiarowy można jeszcze urealnić... I skomplikować. Można zauważyć, że dla 
wielu organizmów, a także dla zjawisk społecznych utrzymanie „status quo” może być 
łatwiejsze niż jego zmiana. Na przykład staramy się postawić namiot w sąsiedztwie dwu 
innych, ale gdy budząc się późnym rankiem widzimy, że pozostał nam tylko jeden sąsiad, lub 
przeciwnie, przybył nam jeszcze trzeci, to machamy ręką i zamiast szukać lepszego miejsca 
idziemy jednak na plażę. Optymalna liczba sąsiadów nie musi zresztą wynosić dwa. Niektóre 
gatunki ptaków do odchowania piskląt potrzebują tzw. pomocników — młodych ptaków, 
często będących rodzeństwem lub wcześniejszym potomstwem gniazdującej pary. Bez 
przynajmniej jednego pomocnika odchowanie piskląt staje się bardzo mało prawdopodobne, 
jednak nadmiar pomocników może stać się kłopotliwy — zaczynają „spiskować” na boku i 
zakładać w pobliżu własne gniazda. Ssak społeczny, np. świstak czy piesek preriowy przeżyje 
przez jakiś czas nawet samotnie, ale może rozmnażać się z sukcesem tylko wtedy, gdy poza 
partnerką ma też kilku dorosłych sąsiadów w kolonii, bo wtedy czas „trzymania warty” przez 
jednego osobnika nie jest za długi i pozostaje mu wystarczająco dużo czasu na poszukiwanie 
pokarmu. Gdy jednak kolonia staje się zbyt liczna narastają konflikty i stresy, a może też po 
prostu zacząć brakować w okolicy pokarmu. Nie koniecznie spowoduje to natychmiastową 


śmierć lub migrację, ale kto rozsądny decyduje się w takich warunkach na posiadanie dzieci? 
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Rysunek 7: Automat komórkowy, w którym komórka rodzi się, gdy ma dokładnie trzech żywych 
sąsiadów, a przeżywa, gdy ma 2, 3 lub 4 żywych sąsiadów. Wycinek po losowej inicjalizacji z 
prawdopodobieństwem 5% po 70 krokach (a) i po 300 krokach (b). 

Tłumacząc to na język reguł automatu komórkowego możemy ustalić, że komórka pozostaje 
„żywa”, gdy ma 2, 3 lub 4 sąsiadów, ale rodzi się — czyli zmienia stan z 0 na 1 tylko wtedy 
gdy ma dokładnie 3 żywych sąsiadów. W takich warunkach samotne pary giną, małe kolonie 
stabilizują się szybko w postaci konfiguracji geometrycznych, w których każda istniejąca 
komórka ma właściwą liczbę sąsiadów, ale w sąsiedztwie nie ma żadnej komórki, która może 
zostać zasiedlona (w terminologii automatów komórkowych struktury o takich 
właściwościach nazywane są „martwymi naturami ). Dopiero odpowiednio liczne kolonie 
mogą dać początek intensywnemu zasiedleniu, które ostatecznie samoorganizuje się w postać 


losowego labiryntu pokrywającego całą dostępną przestrzeń (rys. 7). 


Taki proces wzrostu populacji, który jest początkowo powolny, potem ostro przyśpiesza, by 
na koniec zwolnić, gdy wypełnienie środowiska zbliża się do maksymalnego, powinien 
wydawać się znajomy (rys. 8a). Oczywiście jest to wzrost modelowany ciągłym równaniem 
logistycznym. Zatem model komórkowy okazał się zgodny z matematycznym modelem 
„makroskopowym”, wróciliśmy niemal do punktu wyjścia. Dlaczego ,,niemal” - okaże się gdy 


powiększymy niektóre fragmenty wykresów dynamiki. 
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Rysunek 8: Dynamika populacji zasiedlających obszar 860 tys. komórek wg. różnych modeli 
komórkowych. 23/23 oznacza regułę „żyj gdy masz 2 lub 3 żywych sąsiadów”, 234/234 to reguła „żyj gdy 
masz 2,3 lub 4 żywych sąsiadów” a 234/3 reguła „przeżyj gdy masz 2,3 lub 4 żywych sąsiadów, narodź się 
gdy masz dokładnie 3”. A) Całość historii aż do uzyskania stanu równowagi. B) Przeskalowany końcowy 
odcinek poziomy dla reguły 234/3. C) Przeskalowany końcowy odcinek poziomy dla reguły 23/23. d) Kilka 
pierwszych kroków z uwidaczniające początkowy kryzys liczebności populacji wywołany jej nienaturalną 
strukturą po losowaniu. 

Lokalna, mikroskopowa dynamika nie uśrednia się całkowicie, jak zakłada model 
makroskopowy. Pozostawia swój ślad w postaci minimalnych wahań liczebności populacji w 
stanie równowagi (rys. 8bc), a to, co dzieje się w pierwszych krokach po losowaniu może być 
wyjaśnione tylko za pomocą modelu komórkowego (rys. 8d). W modelach zasiedlania 
wymagających dwu lub więcej „żywych” sąsiadów do wygenerowania nowej „żywej” 
komórki wylosowanie jedynek z pewną gęstością daje populacje o „nienaturalnym” 
rozłożeniu, co powoduje znaczne straty w pierwszym kroku. Dla modeli generujących 
potomstwo w parach straty są łatwo wyrównywane, ale gdy do wygenerowania potomka 
potrzebna jest trójka sąsiadujących zasiedlonych komórek, wzrost populacji rozpoczyna się 


dopiero od czwartego kroku i długo jeszcze pozostaje bardzo umiarkowany”. 
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Drobna modyfikacja, duża różnica 


Na zakończenie musimy zająć się jeszcze jednym modelem, który choć formalnie należy do 
tej samej rodziny, co wszystkie dotąd opisane, jest o wiele bardziej sławny, a jego znaczenie 
wykracza daleko poza modelowanie wzrostu populacji. Można zaryzykować twierdzenie, że 
popularnością wśród laików „bije na głowę” wszystkie inne automaty komórkowe, a nawet 
wszystkie modele matematyczne razem wzięte... Jeśli z modelu opisanego regułą 234/3 
usuniemy możliwość ,,przezywania” komórki otoczonej czterema sąsiadami to otrzymamy 
tzw. „Grę w życie” czyli „Game of Life” wymyśloną przez Johna Hortona Conwaya. Została 
ona opisania przez Martina Gardnera w Scientific American w roku 1970, i od tego czasu 
wciąż fascynuje tysiące naukowców i dziesiątki jeśli nie setki tysięcy hobbystów na całym 
świecie”. 

Dlaczego? „Life” Conwaya nie sprawdza się jako model wzrostu populacji, jest jednak 
znakomitym przykładem układu o właściwościach emergentnych. Z dwu niezwykle prostych 
reguł lokalnych „pozostaw komórkę żywą, jeśli ma 2 lub 3 żywych sąsiadów” oraz „zmień 
komórkę w żywą, jeśli ma dokładnie 3 żywych sąsiadów” układ ten generuje niesamowitą 
złożoność zjawisk. Lokalnie jest on zdolny do osiągania stanu równowagi w postaci zbioru 
różnych, mniej lub bardziej prawdopodobnych „martwych natur” i oscylatorów, ale 
możliwość spontanicznego pojawiania się konfiguracji takich jak „szybowiec? czy LWSS 
(rys. 8), które oscylując przemieszczają się i przenoszą aktywność w obszary już 
ustabilizowane powoduje, że w nieskończonej przestrzeni „Life” nigdy nie mógłby osiągnąć 


stanu równowagi. 


* Internet obfituje w programy implementujące ten model. Chyba każdy fizyk lub matematyk początkujący w 
nauce programowania robi swój Life. Niektóre z tych programów są naprawdę znakomite. Osobiście mogę 
polecić Life32 (http://psoup.math.wisc.edu/Life32.html) zrobiony przez Johana Bontes’a, który dzięki sprytnej 
implementacji uniwersum - „przestrzeni życiowej”, może symulować układy o rozmiarach nawet 


1 mln. na 1 mln. komórek! 
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Rysunek 9: Wycinek inicjowanej losowo przestrzeni w „Game of Life”, będący lokalnie w stanie 
równowagi (a). Widoczne są różne konfiguracje typu „martwa natura” oraz najbardziej popularne 
oscylatory: linia trzech komórek tzw. „blinker” czyli „migacz” oraz cztery „migacze” zestawione razem i 
naprzemiennie tworzące krzyż lub pozbawiony rogów kwadrat. Obok przykłady trzech najbardziej 
popularnych struktur dynamicznych: b) „glider” czyli szybowiec; c) „LWSS” (lightweight spaceship) — inny 
przemieszczający się oscylator; d) „działo szybowcowe” — oscylator produkujący nieskończony ciąg 
szybowców. 

Inwencja badaczy pozwoliła zbadać ogromny zbiór różnych konfiguracji Life, których 
spontaniczne powstanie byłoby bardzo mało prawdopodobne lub wręcz niemożliwe (tzw. 
„rajskie ogrody”). Jedną z pierwszych było tzw. „działo szybowcowe” - oscylator generujący 
nieskończony strumień „szybowców ”. Dzięki niej w listopadzie 1970 r. dowodzony przez 
Billa Gospera zespół z Massachusetts Institute of Technology zdobył ufundowaną przez 
Conwaya nagrodę honorową wysokości 50$ za odnalezienie układu komórek, którego wzrost 
może trwać w nieskończoność. Co więcej odpowiednio skonfigurowane zderzenie kilku 
szybowców może takie działo szybowcowe odtworzyć — czyli „działa są układami 


potencjalnie zdolnymi do samoodtwarzania. Odnaleziono też wiele konfiguracji wędrownych 
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(tzw. „spaceships”) większych i bardziej skomplikowanych niż szybowce oraz konfiguracji, 
które wędrując w nieskończoność zostawiają za sobą Ślad w postaci powtarzalnego ciągu 
innych konfiguracji (tzw. „pufjers”). Wreszcie wynaleziono konfiguracje zdolne do 
wykonywania różnych obliczeń i funkcji logicznych, co nakierowało badaczy na trop 
najdziwniejszej właściwości tego automatu. W roku 1982 udowodniono, że „Gra w życie” 
jest systemem równoważnym uniwersalnej maszynie Turinga, a zatem każdemu możliwemu 
do pomyślenia komputerowi realizującemu algorytmy. Jeśli przyjmiemy, że wszystko na 
świecie da się symulować z dowolną dokładnością to odpowiednio duży „Life” mógłby być 
modelem całego wszechświata. Prawdę powiedziawszy plansza taka musiałaby być 
wielokrotnie większa od symulowanego świata, bo nawet symulacja zachowania 
pojedynczego elektronu wymaga skomplikowanych obliczeń, a nawet proste obliczenia w 
„Life” wymagają tysięcy komórek. Tym niemniej czysto teoretycznie jest to możliwe... Być 
może taki automat byłby nie tylko symulacją, lecz mógłby być nawet równoważny całemu 


światu”, ale tą kwestię musimy pozostawić filozofom i przyszłym pokoleniom naukowców. 


Do poczytania 


e Gardner M. (1970). Games of Life. Mathematical Games in Scientific American 
223(4), 120-123 


e  Sipper M., Reggia A.J. (2001). Idźcie i mnóżcie się. Świat Nauki 10, 54-59 


e Wolfram S. (2002). A New Kind of Science. Wolfram Media, Inc. 


5 Jak na razie nie jest wyjaśnione, czy zjawiska kwantowe dają się modelować algorytmicznie z dowolną 


dokładnością, co stoi na przeszkodzie stwierdzeniu czy wszechświat jako całość jest czy nie jest algorytmiczny. 
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NOTKA: Autor jest z wykształcenia biologiem, przez 10 lat był wykładowcą i asystentem w Instytucie Botaniki 
Uniwersytetu Warszawskiego gdzie zajmował się zastosowaniami w systematyce nieklasycznych metod 
taksonomii numerycznej i analizy obrazu. Jednocześnie (od 1994 r.) w Ośrodku Badania Układów Złożonych 
Instytutu Studiów Społecznych UW zajmuje się modelami  symulacyjnymi procesów 
społecznych, dynamicznych sieci i ewolucji. Po obronie doktoratu w 2001 r. został etatowym pracownikiem ISS, 
a od 2000 r. jest też wykładowcą w Szkole Wyższej Psychologii Społecznej, w szczególności na unikalnym 
interdyscyplinarnym kierunku studiów - Społecznej Psychologii Informatyki i Komunikacji (SPIK). Prowadzi 


różnorodne zajęcia z symulacji komputerowych i seminaria z szeroko pojmowanego ewolucjonizmu. 
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